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ABSTRACT 

In a previous paper, we presented the global solution of a new accretion 
flow model, namely luminous hot accretion flows (LHAFs). In this Letter, we 
first show the corresponding thermal equilibrium curve of LHAFs in the mass 



^ ■ accretion rate vs. surface density diagram. Then we examine its thermal stability 

again local perturbations. We find that LHAFs are thermally unstable when 
thermal conduction is neglected. However, when the accretion rate is not very 
large, the timescale of the growth of perturbations is longer than the accretion 
timescale, therefore the instability has no dynamical effect on the accretion flow. 
When the accretion rate is large, the perturbations can grow very fast at a certain 
radius. As a result, some cold clumps may form and the accretion flow will 
q . become multi-phase. 

CO 

Subject headings: accretion, accretion disks — black hole physics — galaxies: 



active — galaxies: nuclei 



1. Introduction: the physics of the new hot accretion disk solution 

There has been a great interest to the accretion process around black holes. The most 
famous accretion solution is the geometrically thin and optically thick cold disk model de- 
veloped by Shakura & Sunyaev (1973; hereafter SSD) and others. The second solution was 
discovered by Shapiro, Lightman & Eardley (1976; hereafter SLE). This solution is optically 
thin and hot, with T e ~ 10 9 K. Pringle (1976) found that SLE is thermally unstable although 
it is not clear what is the consequence of such instability. 



In both SSD and SLE, the energy advection is neglected. It was found that when the 
mass accretion rate is higher than the Eddington rate M Edd (= 10L Edd /c 2 ), the large optical 
depth traps most of the photons therefore most of the viscously dissipated energy is stored in 
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the gas and advected into the black hole rather than radiated away ( "slim disk" ; Abramowicz 
et al. 1988). 

The fourth accretion solution is the advection-dominated accretion flow (ADAF; Narayan 
& Yi 1995; Abramowicz et al. 1995, hereafter A95; see reviews by Narayan, Mahadevan & 
Quataert 1998; Kato, Fukue & Mineshige 1998). Different from SLE, the energy advection 
is included in the ions energy equation of an ADAF, Q adv = Qvis — Qie- Here Q adv , Qvis and 
Q lc are the rates of energy advection, viscous heating and Coulomb cooling per unit area 
of the accretion disk, respectively. In a typical ADAF, the density of gas is very low, so 
Q\c <C Qv\s ~ Qadv, i-e., the viscous heating is balanced by advective cooling. Since Q 1C oc M 2 
while Q v i S oc M, i.e., Q ie increases faster than Q vis with increasing M, there exists a critical 
rate Mi, determined by Q vis ~ Q ic . At this rate, a large fraction of the viscously dissipated 
energy is transferred to the electrons and radiated away, so the accretion flow ceases to be 
advection-dominated. Above Mi, it was thought previously that no hot accretion solution 
exists and the only viable solution is SSD. 

However, our recent work (Yuan 2001, hereafter Y01) indicated that this is not true: 
above Mi, a new hot accretion solution exists. To illustrate, let's first write out the formula 
of the energy advection of ions, 

k 1 dT kT I dp' 



Qadv = PHV T 



Qint Qcom; (1) 



pm^ 7 — 1 dr pm^ p dr 

i.e., the energy advection consists of two terms, namely the internal energy gradient term 
Qint and the compression work Q com - The full condition for the existence of a hot accretion 
solution, when the flow starts out hot, is dT/dr < 0, i.e., Q int > (note v r < 0). For an 
ADAF, M < Mi, Q adv = Q vis - Q ie > 0, so Q int = Q adv + Q com > 0. When M > M, Q adv = 
Qvis — Qie < 0, but there obviously exists another critical accretion rate, M2, determined by 
Qcom + Qvis ~ Qie- Below M 2 , we still have Q int = Q vis + Q CO m - Qie > 0. Therefore, the 
accretion flow can remain hot if it starts out hot. We denote this new solution as luminous 
hot accretion flow (LHAF hereafter). 

From the above analysis, we see that LHAFs are along the line of ADAFs — the equations 
describing both are completely the same and we just extend ADAFs to higher accretion rates. 
However, the dynamics of an LHAF is quite different from an ADAF. In an ADAF, Q adv > 0, 
i.e., the advection plays a "cooling" role from the Lagrangian point of view. But in an LHAF, 
Qadv < 0, so it plays a "heating" role. In the language of entropy, it is the conversion of 
entropy together with viscous dissipation that supplies the radiation of LHAFs. In this sense, 
an LHAF is dynamically similar with Bondi accretion and cooling flow in galactic clusters. 

In this context we can understand why previous authors didn't find this solution. In 
Narayan & Yi (1995) and Esin et al. (1997), they a priori set the advection factor "/"(= 
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Qadv/Qvis) as positive, so they only obtained the ADAF solution. Similarly, A95 and Chen 
et al. (1995, hereafter C95) didn't find this solution because they set the parameter £ (see 
eq. (6) below) as positive. 

In this Letter, we will first show the thermal equilibria of LHAFs in the M vs. £ (surface 
density) diagram (§2). Such a diagram is widely used in the study of accretion disks (e.g., 
A95; C95; Kusunose & Mineshige 1996; Bjonsson et al. 1996). Then in §3 we investigate 
the thermal stability of LHAFs. 



2. Thermal equilibrium curve of LHAFs 

We take a one-temperature accretion flow as an example for simplicity. For the purpose 
of comparison, we adopt almost exactly the same equations as C95 (see also A95). 



M = -2irREv r , (2) 

= |W\ (3) 

QaAw Qvis Qr&A- (4) 

They represent conservations of mass, angular momentum, and energy, respectively. Here 
v = ^ac s H is the kinetic viscosity coefficient, E = 2pH is the surface density, /* = 1 — 

9tt(3R s )/ltt(R)(R/R s ) 2 ] with R s = 2GM/c 2 ,g = -f (dlnQ/dlnR). The forms of Q vis , Q adv 
and Q rad are, 

Qvis = ——Q 2 f*g, (5) 



47T 

p 2irR 2 p"- 



Q adv = -Sw r -£ = (6) 



Q rad = 8 ( xT 4 f^ + v / 3 + ^y . (7) 

\ ^ Vbrem / 

The optical depth r = r cs +r abs , with r cs = (l/2)/c cs S and r abs = Qbrem/ (8cxT 4 ). The equation 
of state is p = p gas + p rad , with p gas = ^pT, and p rad = + -^j (Abramowicz et al. 

1996). The bremsstrahlung radiation Qbrcm = ^H(q ei + q ee ). The formula for q ei and q ee are 
from Narayan & Yi (1995). Another useful equation is the hydrostatic balance equation, 
H = c s /VL K , with c s = (p/p) 1 / 2 is the local sound speed. The Paczyhski & Wiita (1980) 
potential is used. 
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We solve eqs. (2)-(7) to obtain the thermal equilibrium curve of LHAFs. One noticeable 
parameter is £ in eq. (6). From §1, we know that for an LHAF, Q a dv < 0, therefore we 
should set £ < to recover this solution. We can recover all other accretion solutions by 
simply setting £ = 1, as in A95. The crudeness of the value of £ doesn't affect our qualitative 
results although it does prevent us from obtaining the exact quantitative results such as the 
ranges of accretion rates to which various accretion disk models correspond. 

Figure 1 shows the thermal equilibrium curves of different solutions. Comparing with 
the figures in A95 or C95, we see that the three "U" shaped thick lines are new and they 
denote LHAFs with different values of £. These lines extend from optically thin regime to 
optically thick regime, bridging the SLE and the SSD solutions. It is not clear what physical 
implications the "bridging" has, or it might just be a mathematical trick. But we note in this 
context that in the slim disk solutions, they did find solutions with the advection term being 
negative in some cases (Abramowicz et al. 1988; Szuszkiewicz, private communication). 
Some segment of the "U" curves are superimposed on the SLE and SSD lines. We find 
that advection in these "superimposed" segment is equal to zero, same with SLE and SSD. 
Only the "separate" segment denotes the genuine LHAFs, since for this segment the energy 
advection is negative and plays a significant role in the energy balance of ions. We see that 
above the critical accretion rate of ADAF, Mi(~ 0.1M Edd ), both ADAF and SLE disappear 
and LHAF is the only available hot solution. This is consistent with our analysis in §1. 

For comparison, we also solve the global solutions of the standard accretion equations 
describing a one-temperature accretion flow with the same parameters as in Fig. 1 (see 
Yuan 1999 for the equations except that the radiative term is now replaced by eq. (7) of 
the present paper). The results are denoted by the filled circles (for ADAFs) and triangles 
(for LHAFs). We can see that our local algebraic analysis is qualitatively consistent with 
the global results, i.e., when the accretion rate is lower than a certain value, the solutions 
are ADAFs; while above this value, the solutions are LHAFs. We also calculated the value 
of £ for a specific example of M — 10~ L6 MEdd and the outer boundary conditions of R ont = 
10 3 R s ,T out = 2 x 10 9 K, and v/c s = 0.4. We obtained £ « -0.3 &tR = 5R S . 

We want to emphasize that our "U" shaped line is completely different from the line 
with the similar shape in C95, which was obtained by setting a very large viscous parameter 
a > a cr with a cr ~ 0.2 or larger. In fact, a later more accurate treatment of microphysics and 
the inner boundary condition by Bjornsson et al. (1996) indicated that a cr > 1, therefore, 
the "U" shaped branch in C95 is unphysical. 
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3. Thermal stability 

From the density profile of the global solution in Y01, we know that LHAFs are viscously 
stable. An important problem then is to analyze its thermal stability. Before doing that, we 
note that an LHAF is dynamically very similar to a cooling flow in galactic clusters. Many 
authors have studied the thermal stability of cooling flow and concluded that it is unstable 
against local perturbations if thermal conduction is neglected (e.g., Fabian & Nulsen 1977; 
Mathews & Bregman 1978; Nulsen 1986). Thermal conduction can strongly stabilize the 
cooling flow (Zakamska & Narayan 2003; Kim & Narayan 2003). 

Now let's investigate the thermal stability of LHAFs against local perturbations. Such 
perturbations can be, e.g., that the local gas density increases wherever the initial ambient 
magnetic field is less than average. We follow the standard analysis approach presented by 
Kato, Abramowicz & Chen (1996). Taking the surface density as an example, we denote 
perturbations Si superimposed over the unperturbed quantities S as a = Si/S oc exp(nflt — 
ikr) where nVL is the growth rate of perturbation and k is the perturbation wave number. 
Substituting the perturbed quantities to the time-dependent accretion equations, we can 
obtain the perturbed equations and the dispersion relation. The arguments for ADAFs 
presented in Kato et al. (1996) hold for LHAFs as well although the advection is negative 
and radiative cooling is important here. We assume that thermal conductivity is suppressed 
by the presence of the tangled magnetic field in the accretion flow (but see Narayan & 
Medvedev 2001). We expect the flow will be thermally stable if thermal conduction is 
strong. The dispersion relation is, 

-3ntta = Ga, (8) 

and the condition for instability is 

G _ Qvis ( d\nrj \ Qrad ( dlllQrad 

W V<9hiSy w W V <9hrS 

Here W = 2Hp and r\ = z/S. Since r\ oc Sc 2 /f2 oc W/Q, so (dlnrj/dlntyw = 0. For a 
one-temperature LHAF, Q rad oc Hp 2 T l l 2 oc Y?T l l 2 jE oc S 2 , so 

(<91nQ rad /,91nSV > 0. (10) 

Therefore one-temperature LHAFs are thermally unstable against local perturbations. 

The case of two-temperature LHAF is a little bit more complicated. Treating the energy 
equation in Kato et al. (1996) as the energy equation of ions, the instability condition is 
exactly the same as eq. (9) except that we should replace Qrad with Q ie , the Coulomb 
collision cooling rate. Since Q ie oc p 2 ifT e ~ 3 ^ 2 Tj, to determine the sign of G, we must turn to 
the electron energy equation to obtain the relation between T e and T { . Writing the energy 
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equation of electrons as Q a dv,e = Qie — <5rad,c, with Q a dv,e and Q ra d,e are the energy advection 
and radiative cooling of electrons respectively, the global solution of Y01 indicates that for 
LHAF a good approximation to the electrons energy equation is Qi e ~ Q a dv,e hi the region 
outside of ~ 100R S and Qi e ~ Q r &d,e inside 100R S (see Figure 5 in Y01). Taking Q ie m <5adv,c 
we can get T e oc W 1 ^ 1 ^. Hence Q ie oc E 6 / 5 ^ 1 / 5 and 



So a two-temperature LHAF is thermally unstable at least at the region outside of ~ 100i? s . 
Inside ~ 100R S , the radiative cooling is in general dominated by the thermal Comptonization. 
The seed photons can be the synchrotron photons, or other soft photons from cold matter 
such as a SSD underlying or outside of the LHAF, or cold clumps suspending in the hot gas. 
The origin of the cold clumps in LHAFs can be due to the thermal instability of LHAFs (see 
below), or originally existed in the accretion material. For luminous sources the soft photons 
from the cold matter is likely to dominate over the synchrotron photons as the seed photons 
of Comptonization (e.g., Zdziarski et al. 1998). Since Q ra ,d,c = (4kT e / m e c 2 ) pH K es U r c, with 
U r denotes the flux of seed photon, using Q ie Q ra ,d, e , we find that if U r oc Tf or stepper, the 
two-temperature LHAFs will be thermally stable, otherwise it is unstable. Unfortunately, 
it is not easy to determine whether this condition can be satisfied or not, because of the 
complicated energy interaction between the cold and hot phases in accretion flow (e.g., 
Ferland & Rees 1988). 

Now let's discuss the consequence of the thermal instability against local perturbations. 
An important factor is how fast the perturbations grow compared to other related timescales. 
It is believed that as the result of instability cold dense clumps will form in the hot flow 
if the growth of perturbations is fast enough. This mechanism has been used to explain 
the condensations of galaxies from the intergalactic medium, formation of solar prominence, 
and condensations in planetary nebulae (e.g., Field 1965). In the case of accretion flows, we 
should compare the accretion timescale t acc with the growth timescale of perturbations t grow . 
For LHAFs, from eq. (8), the timescale of the growth of local perturbations is 




(11) 




(12) 



This equals the thermal timescale of accretion flows, t t h, 



t. 



3(7-1) 



(13) 



grow 



2 



for adiabatic index 7 = 5/3. 
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We compare the two timescales, t acc and t gr0 w, by numerical calculations. This requires 
us to solve the global solution of two-temperature LHAFs. This was done in Y01. The 
results indicate that depending on the value of M, there are two types of LHAF. When 
Mi < M < M 2 , with M 2 can be as high as 5 Mi, the accretion gas in LHAF is hot throughout 
the disk; when M > M 2 , up to Mgdd; within a certain radius the density of the accretion gas 
is so high that even the sum of compression work and viscous dissipation can not balance the 
strong Coulomb cooling. As a result, the hot accretion flow will collapse onto the equatorial 
plane and form an optically thick cold annulus. While the quantitative results such as the 
value of M 2 depend on parameters such as a, the qualitative results should not. 

Correspondingly, there are two kinds of results when we compare t acc and t gr0 w, as shown 
by the two lines in Figure 2. The corresponding two models are taken from Fig. 4 of Y01. 
The solid line corresponds to M = 0.1M E dd < M 2 , so the accretion flow is hot throughout the 
disk; the dashed line to M — 0.3M E dd > M 2 , so the hot accretion flow collapses at a radius 
and form a cold annulus. We see that when M < M 2 , t S[OW /t a/CC > 1, so the perturbations 
have no time to grow before they are swallowed by the black hole, therefore, the thermal 
instability will have no effect on the dynamics of LHAFs. When M > M 2 , t grow /t acc < 1 
at the transition radius. As a result, some cold dense clumps will form there if LHAFs 
are thermally unstable at that region. We speculate that once the condensation of the hot 
phase makes t gr0 w ~ ^acc, the condensation will stop. This constraints will be useful when 
quantitatively calculating the physical states of clumps such as their filling factor, density 
and temperature. 

In our calculations of t acc and t gr0 w, we considered only the "standard" radiative pro- 
cesses, namely synchrotron, bremsstrahlung, and their Comptonization. If we include the 
additional Compton cooling due to the soft photons from the cold phase, t grow will decrease ( 
Mi and M 2 will also decrease). We expect that we will still have t grow > t acc when M < M 2 . 

Note that LHAF will not collapse into an optically thick cold disk as the result of local 
perturbation since the perturbation wavelength is smaller than the disk scale- height. The 
thermal stability of LHAFs against long wavelength perturbations is hard to determine, 
and it is also not clear how fast the growth of perturbations is compared to the accretion 
timescale if the flow is unstable. On the other hand, if the filling factor of the cold clumps 
is large enough, it is possible that these clumps may assemble due to frequent collisions and 
form a disk-like large scale structure. 

Many authors have proposed such cold/hot two-phase accretion flow model for AGNs 
(e.g., Guilbert & Rees 1988; Ferland & Rees 1988; Kuncic, Celled & Rees 1997). But almost 
all these work focus on the thermal state of cold clumps and how the clumps re-radiate the 
energy they absorb. Krolik (1998) also suggested a two-phase accretion flow model. Different 
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from our suggested physical mechanism of clumps formation, he suggested that such a two- 
phase feature is the result of the instability in the radiation-pressure-dominated innermost 
region of a SSD (see e.g., Gammie 1998). 

I thank Shoji Kato, Andrei Beloborodov, and Ramesh Narayan for valuable discussions, 
and the anonymous referee for helpful comments. This work was supported in part by NASA 
grant NAG5-10780 and NSF grant AST 0307433. 
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Fig. 1. — The thermal equilibrium curve of various accretion solutions. The accretion rate 
is in units of M-^dd = 10-^Edd/c 2 and the units of E is g cm -2 . The parameters are M/M & = 
10, a = 0.1, and R = 5R S . The thin solid lines are for £ = 1, representing ADAF, SLE, 
SSD, and slim disk. The thick solid, dotted and dashed lines are for £ = —1, —0.1 and —10, 
respectively, representing LHAFs. The squares are the location of unity scattering optical 
depth. The filled circles and triangles are the results obtained from the global solutions of 
ADAFs and LHAFs, respectively. 
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Fig. 2. — The ratio of the accretion timescale to the timescale of growth of perturbations as 
a function of radii for two LHAF solutions presented in Fig. 4 of Y01. The solid (dashed) 
line is for M = 0.1(0.3)M Edd . Other parameters are a = 0.1,(3 = 0.5. 



